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SUMMARY 


The profiles and thicknesses of normal shock waves in argon at Mach 
numbers of 1-335, 1.454, l-5?6, and 1-713 were determined experimentally 
by means of a free-molecule probe whose equilibrium temperature is related 
by kinetic theory to the local flow properties and their gradients. 

Comparisons were made between the experimental shock profiles and 
the theoretical profiles calculated from the Navier-Stokes equations, the 
Grad 13-moment equations, and the Burnett equations. New, very accurate 
numerical integrations of the Burnett equations were obtained for this 
purpose with results quite different from those found by Zoller, to whom 
the solution of this problem is frequently attributed. 

The experimental shock profiles were predicted with approximately 
equal success by the Navier-Stokes and Burnett theories, while the 
13-moment method was definitely less satisfactory. A surprising feature 
of the theoretical results is the relatively small difference in predic- 
tions between the Navier-Stokes and Burnett theories in the present range 
of shock strengths and the contrastingly large difference between predic- 
tions of Burnett and the 13-moment theories. It is concluded that the 
Navier-Stokes equations are correct for weak, shocks and that within the 
present shock strength range the Burnett equations make no improvement 
which merits the trouble of solving them. For shocks of noticeably 
greater strength, say with a shock Mach number of more than 2.5, it 
remains fundamentally doubtful that any of these theories can be correct. 


INTRODUCTION 


A new experimental method for measuring shock-wave structure involving 
the use of a free-molecule probe was introduced by Sherman in reference 1. 
He made his measurements in a low-density wind tunnel, where the shock 
waves produced were thick enough to permit local measurements within them 
while the probe, being small compared with the mean free path, did not 
disturb the macroscopic flow pattern. This technique enabled him to 
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obtain information on the shape and thickness of the shock profiles which 
could be compared in a detailed fashion with the predictions of several 
theories for shock structure. 

The free -molecule probe used by Sherman was essentially a cylindrical 
wire whose temperature could be measured. The wire was stretched across 
the jet of the wind tunnel and was mounted on a mechanism which traversed 
it through the normal shock waves formed at the open end of a hollow con- 
ical shock holder. The data obtained consisted of point -by-point profiles 
giving the variation of wire equilibrium temperature with position through 
the shock waves. 

Experimental wire-temperature profiles were compered with predictions 
based on the Navier-Stokes equations, the Grad 13-moment equations, and a 
kinetic-theory method of Mott-Smith. The Burneot equations were not used 
by Sherman because of the difficulty in solving them for the shock-wave 
boundary -value problem. 

Sherman ! s experimental data were much better predicted by the Navier- 
Stokes equations than by the more complex kinet:.c-theory methods. However, 
the comparison among the various theories was somewhat clouded by two 
facts: (l) Most of the experimental data was obtained for diatomic gases, 

whereas the kinetic theory has been fully developed only for monatomic 
gases, and (2) the shocks tested were too strong, falling outside the range 
of Mach numbers for which a single-valued shock-transition solution of the 
13 -moment equations exists. 

In an attempt to remove the two objections listed above, the meas- 
uring technique of reference 1 was applied to shock waves in argon at 
Mach numbers of 1.335* 1-^5^* 1*576, and 1.713* On the basis of the 
experimental accuracy obtained by Sherman and t3ie magnitude and character 
of the differences in predictions between the Navier-Stokes and 13-moment 
theories, it was hoped that these new experiments would be capable of 
discriminating clearly between the competing theories. The present report 
describes the new measurements and rounds out the comparison with theory 
by including high-speed computing-machine solutions of the Burnett 
equations . 

This investigation was conducted at the University of California 
under the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. The authors are also grateful to 
the Naval Research Laboratory, where the numerical integrations of the 
Burnett equations were performed. 
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SYMBOLS 

Zoller*s shock-strength parameter 
coefficients in equation (A5) 
dimensionless peculiar molecular velocity 
specific internal energy, cm^/sec^ 
molecular velocity-distribution function 
Bell-Schaaf speed-ratio functions 

convective heat-transfer coefficient for free- 
molecule flow 

modified Bessel function of first kind, order zero, 
and argument s^/2 

modified Bessel function of first kind, order one, 
and argument s^/2 

electric current, amp 

indexes, used as subscript (i,j = 1,2,3) 
thermal conductivity of wire, ergs/cm sec °K 
thermal conductivity of gas, ergs/cm sec °K 
reference length, p*ypu 

reference length, p-^pu (used in table I) 
half length of wire 
shock Mach number 

mass flow constant, equations (B3), grn/cm^ sec 
number density of gas molecules, cm“3 
momentum constant, equations (B3)j dynes/cm^ 
thermodynamic pressure, dynes/cm^ 
impact pressure, dynes/cm^ 
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components of viscous-stress tensor, e.g., Pxx^ 
Pyy; P zz > dynes/cm 2 

components of excess viscous-pressure tensor, 
dyne s / cm2 

Navier-Stokes part of p , dynes/cm^ 

Burnett part of p.^, dynes/cm^ 
reservoir stagnation pressure, dynes/cm^ 

energy constant, equations (B3), ergs/cm^ sec 
heat-flux vector component, e.g., q x , ergs/cm^ sec 

heat flux, ergs/cm^ sec 
Navier-Stokes part of q y 
Burnett part of q x 

gas constant for argon, ergs/cm °K 
resistance of wire per ur.it length, ohms/cm 
resistance of wire at terperature T 0 , ohms 

resistance of wire at temperature T w , ohms 
functions defined in equations (Bll) 
radius of wire, cm 
molecular speed ratio 
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gas temperature, °K 

adiabatic wire equilibrium temperature, °K 
stagnation (room) temperature, °K 
wire equilibrium temperature 

2 

reference temperature, °K, T 0 

7+1 

gas velocity, cm/sec 
mean flow velocity, cm/sec 
dimensionless velocity, equations (B8) 
boundary value of v 

normalized, dimensionless velocity, equations (Bl8) 
distance along flow, cm 

Grad's dimensionless distance, equation (l) 

Zoller's dimensionless distance, equations (B8) 
thermal accommodation coefficient of wire 
temperature coefficient of resistance, °K“^* 
specif ic-heats ratio 

maximum-slope shock-wave thickness, based on veloc- 
ity profile and measured in units of x and y, 
respectively 

maximum-slope shock-wave thickness based on wire- 
temperature profile 

Grad's shock-strength parameter, equation (2) 
radiant emissivity of wire at temperature T 0 



radiant emissivity of wire at temperature T v 

dimensionless Navier-Stoi.es stress, equations (B8) 

dimensionless Burnett stress, equations (B8) 

initial stress perturbation amplitude, equations (Bl4) 

temperature perturbation, equations (B13) 

initial temperature perturbation amplitude, equa- 
tions (bi4) 

Burnett heat-flux coefficients 

Maxwellian mean free path ahead of shock, cm 

root of characteristic equation (Bl6) 
viscosity, poises 

viscosity at reference temperature T*, poises 
viscosity ahead of shock; poises 

boundary value index, equations (B12) 
absolute molecular velocity, cm/sec 
ith component of | 

dimensionless Navier-Sto tes heat flux, equations (B8) 
dimensionless Burnett heat flux, equations (B8) 

initial heat-flux perturbation amplitude, 
equations (Bl4) 

gas density, gm/cm^ 

Stephan -Boltzmann constait, ergs/cm^ sec °K 
dimensionless temperature variable, equations (B8) 
boundary value of t 

velocity perturbation, equations (Blj) 
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fi initial velocity perturbation amplitude, equa- 

tions (Bl4) 

cjo exponent of viscosity-temperature lav 

on Burnett stress coefficients, used in appendix B 

EXPERIMENTAL APPARATUS 


Mach 1.5 Nozzle 

The experiments were performed in the open-jet, continuous -flow, 
lov-density wind tunnel at the University of California, Berkeley. This 
wind tunnel is described in reference 2. For the purposes of the present 
experiments a nozzle was designed and built to give an exit Mach number 
of 1.5 in a monatomic gas at design flow conditions. This proved to be 
an extraordinary supersonic nozzle. Evidently, because the combination 
of low exit Mach number and extremely thick nozzle boundary layers, the 
nozzle is not subject to the formation of oblique shocks or sharp rare- 
faction fans at the exit lip when the pressure in the test chamber is 
raised or lowered from the exit static, pressure of the open jet. What 
appears to happen instead is that a rise in test-chamber pressure at any 
given stagnation-chamber pressure only produces a smooth increase in 
boundary -layer thickness along the entire nozzle length. This in turn 
yields an effective nozzle of smaller bore and area ratio than the orig- 
inal but of equal contour smoothness. Lowering the chamber pressure just 
reverses the effect. The outcome of this behavior is a great flexibility 
of selection of Mach and Reynolds numbers, obtained by independent adjust- 
ment of stagnation- chamber and test-chamber pressures, without sacrifice 
of flow uniformity. This was exploited in the present experiments by 
making all shock-wave surveys at approximately the same free -stream mean- 
free-path length. The four shocks were then of roughly equal thickness 
and could be surveyed with correspondingly equal accuracy. The nozzle 
exit diameter was 7*91 inches, and the diameter of the isentropic core 
of flow varied from 3 to 4 inches over the range of Mach numbers covered 
here. Lop act pressure surveys taken at the exit plane are shown in 
figure 1. 


Nozzle Calibration 

Because of the extreme importance of an accurate determination of 
Mach number for the present tests, the method of nozzle calibration was 
very critical. Several different independent methods for determining 
the Mach number were used, and the worth of the final calibration is 
shown by the excellent agreement obtained among these methods. 
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The instruments used were a 0.300-inch-dianeter source-shaped impact 
probe, corrected for viscous effects by the methods of reference 3, and 
a 0.300-inch-diameter static -pres sure probe. The latter had 8 pressure 
orifices equally spaced around its periphery at a cross section 13 inches 
from the 10° conical tip and 13 inches upstream from the cross-stream 
support. When the pressure orifices were positioned in the test region, 
the tip of the probe protruded into the stagnation chamber so that tip 
effect was avoided. The 13-inch length downstream of the pressure orifices 
permitted the probe to be traversed throughout :he axial length of the 
nozzle and minimized cross-stream support interferences. The static probe 
was assumed to be free of viscous effects. Its diameter was small com- 
pared with the nozzle throat diameter of approximately 6 inches. No 
change in nozzle wall pressure was observed upoi insertion or removal of 
the probe, and excellent agreement was found between wall pressure and 
probe pressure under conditions where the probe surveys indicated constant 
static pressure across the jet. 

The two local pressure measurements, impact p^ and static p, were 
combined to determine the local Mach number from the ratio P^/P- In 
addition, the Mach number was determined from tie ratios Pi/P 0 and 
p/p , where p Q is the stagnation pressure, assuming an isentropic flow. 

As a final check on the internal consistency of the Mach numbers 
determined in the above fashion, the free-strean equilibrium temperature 
of the free-molecule probe used for the shock-wave surveys was measured. 
Although the probe was initially uncalibrated, in that its radiation 
losses were unknown but significant, the probe «ras extremely sensitive 
to Mach number variation in the Mach number range of the tests. A smooth 
correlation between probe equilibrium temperatures and Mach numbers deter- 
mined from pressure measurements was essential for confidence in the 
latter . 

The final calibrated flow conditions used in the tests are listed 
in table I. The Mach number values, cited to three decimal places, repre- 
sent averages of the values obtained from the several methods. Individ- 
ual values differed from these averages by less than 1 percent. 


Shock Holder 

A single shock holder, similar in design to those employed in the 
tests of reference 1, was used throughout the present experiments. It 
is shown in figure 2. The entry diameter was 3 inches, which is about 
the dimension of the isentropic core of the nozzle flow at the lowest 
Mach number M]_ = 1.335 (see fig. l). The shock holder was slotted to 
provide for motion of the wire probe downstream of the entry plane, but 
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this provision turned, out to he unnecessary since all the shock waves used 
in the tests were slightly detached. An internal included angle of diver- 
gence of 40° was selected for the shock holder, which represented a com- 
promise between the estimated minimum internal angle required to allow for 
boundary -layer growth within the shock holder and the maximum external 
angle for an attached shock wave. These two conditions could not be met 
simultaneously at the low Mach numbers of the tests, and therefore the 
shock waves were slightly detached. However, the glow photographs of the 
shock-holder operation (fig. 3) indicate that at least the central por- 
tions of the shock waves were plane. The shock- detachment distance 
decreased with increasing Mach number. 


Free -Molecule -Flow Equilibrium-Temperature Probe 

The 0.00025 -inch-diameter tungsten resistance-wire probe described 
in reference 1 was used for all the shock-wave measurements. A schematic 
diagram of the resistance-wire instrument is shown in figure 4. In every 
case the Knudsen number of the probe (the ratio of the mean free path in 
the gas to the probe diameter) was greater than 50, which insured that 
free-molecule-flow conditions always prevailed. 

The wire resistance was found to vary linearly with temperature over 
the temperature range encountered in the tests . An oven calibration gave 

(R') w = (R') 0 [l + o^T*. - T 0 )J 

with (R') 0 equal to about JO ohms at 23° C (the exact value differing 
from probe to probe, depending on the length of the measuring section) 
and cap = 0-00332 per °C, also at 23° 0 - 

The wire current of 0.100 ± 0.001 milliampere was supplied by a 
q.^b-volt mercury cell in series with dropping resistances and was meas- 
ured by observation of the potential drop across a 30.0-ohm resistance. 

This potential drop and also that across the measuring section of the 
freei-molecule probe were measured with a Rubicon hand-balancing potenti- 
ometer and a Leeds and Northrup galvanometer. The combined potentiometer - 
galvanometer system had a resolution of better than 1 microvolt (0.1 micro- 
volt could be estimated)^ which is about 0.5 percent of the probe potential 
change between upstream and downstream limits of the weakest shock wave 
tested. Alinement of the free-molecule wire parallel to the plane of the 
shock-holder entry was accomplished by adjusting the wire supporting frame 
and by bending the heavier "tails” of the wire leads. The alinement was 
checked before and after the tests, and it was estimated that the maximum 
misalinement error was not greater than about 0.003 inch over the length 
spanned by the measuring section. If it is assumed that this misalinement 
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only causes the wire to read the average signal, of a perfectly alined wire 
over a 0.0015-inch-travel through the shock, one can demonstrate graphically 
that this produces no significant change in the measured shock profile. 


Wire-Traversing Arrangement 

The wire setup, shown in figure 2, is essentially the same as the 
system described in reference 1. The wire was moved by a small electric 
motor; its position was indicated by a 1-inch-nravel, 0.0001-inch-reading 
Ames dial indicator. 


Flow Visualization 

Unfortunately, an extremely bright afterglow in argon was discovered 
only after the program was completed and, consequently, did not contribute 
to the exploratory phases of the experiments. However, pictures were 
obtained with it (fig. 5) which confirmed, a posteriori, certain important 
conclusions about shock placement and shape wh:.ch had been formed from 
observation of wire -temperature traces. 


Gas Supply 

The argon used for the tests was supplied by Linde Air Products Co. 
and was guaranteed to be 99*9 percent pure. It was introduced directly 
into the wind tunnel without additional drying or purification. Some air 
leaks into the wind tunnel, of course, and may mix with the argon in small 
amounts. What one fears is the presence of an impurity of a molecular 
weight very different from that of the test gar. in sufficient amount to 
add diffusion to the dissipative mechanisms broadening the shock. It can 
be shown theoretically (e.g., ref. 4) that 1 percent of air in argon, a 
very liberal allowance for possible leakage, produces less than l/2 per- 
cent shock thickening by diffusion. 


EXPERIMENTAL RESULTS 


Free-Stream Probe Readings and Radiation Corrections 

Initial measurements were made with the free-molecule probe in the 
uniform stream at the Mach numbers chosen for the final shock-wave test. 
These data are shown in figure 5* The wire -temperature readings fell 
consistently below the theoretical values given by the analysis of ref- 
erence 5* On examining the possible reasons for the discrepancy it was 
found that the losses due to heat conduction were small but that a sizable 
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loss might be expected from radiative heat transfer. An acceptable fit 
to all the free-stream data was obtained with a value of e 0 /a, =0-13 
where e Q is the emissivity of the wire at room temperature and a is 
the thermal accommodation coefficient of argon on tungsten (see 
appendix A) . 

The existing experimental arrangement did not permit accurate inde- 
pendent determinations of the quantities e G and a. However, the cor- 
relation between free-stream Mach number and probe temperature, using 
the single value e 0 /a = 0 . 13 , was sufficiently good to warrant acceptance 
of this value as an accurate calibration constant for the probe. 

It should be noted that the formula for free -molecule-flow convective 
heat transfer given by equation (A2) of reference 1 is in error. The cor- 
rect formula is given by equation (A3) of this paper. 


Shock-Wave Profiles 


The principal results of the experiments are shown by the data points 
plotted in figures 6(a) to 6(d). In these figures the ratio of the meas- 
ured wire temperature T w to the measured stagnation temperature T 0 is 
plotted against the nondimens ional space variable y. This nondimensional 
length is related to the physical distance x of the shock-wave traverse 

by 


y 


it ex 

35 L* 


( 1 ) 


in which e is the shock strength introduced by Grad (ref. 6) and is 
defined in terms of the initial Mach number M-j_ by 



3 + 5M L 2 


and. L* 
cosity 


is a reference length pertaining to the 
at the temperature — T 0 = ^ T 0 . 


value jx* of the vis- 
it is defined by 



( 3 ) 
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The value T* = 


corresponds to the va _ue of the gas temperature 


7+1 ° 

at the sonic point of an adiabatic flow having :otal temperature T Q . 


The major difference between the method used here to present the data 
and that employed in reference 1 is that in this report the data have not 
been "normalized" so as to make the downstream end points of the shock 
wave agree with theory. The normalization technique used in reference 1 
is open to some criticism; although it does not alter the maximum-slope 
thicknesses of the shock waves , it obscures some features of the data. 


THEORETICAL CALCULATIONS 


The theoretical procedures for finding shock-wave profiles in terms 
of the temperature of a free-molecule wire with no radiation or metallic 
conduction losses are given in detail in reference 1 for the Navier-Stokes 
and 13-moment equations. Extensions of these procedures to correct the 
predicted profiles for radiation losses and to Include use of the Burnett 
equations are given in appendixes A and B of th* present report. The fol- 
lowing paragraphs will review, for convenience, only the basic philosophy 
of the theoretical approach. 

The computation of free-molecule-flow energy transfer to the wire 
requires knowledge of the distribution of molecular velocities at the 
wire location. Thus, the initial object of the analysis is to find the 
molecular distribution function f at each point within the shock wave. 
For the calculations discussed here f is assuned in the form of a local 
Maxwellian distribution skewed by the viscous soress and heat flux; that 
is, 


f(ijx) = 


n 


.-C2 


(2«RT) 


3/2 


1 + 


p i; 


CjCh 


q.C. /' 
SLA_; 2 c 2 


p '2RT 5 


\ 

- l! 
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In thin C is the peculiar molecular velocity absolute velocity | 

\ 


minus r 


u, 


:;ean flow velocity u ) divided by \[2R T. All the parameters ( n, 

m _ \ o >-» c -Tiinoti one v ' 


P> T, 


P,- v 


and q.j ^ are functions of x. 


If this distribution function is substitutid into the Boltzmann 
equation and the moments of that equation taken with respect to 1, |j_, 

and - E.f., one obtains the familiar conservation equations of fluid 
2 11 

mechanics . These equations initially form an indeterminate set because 
of the explicit appearance of the stress- and hsat-flux variables. The 
problem at hand is then threefold: (l) The conservation equations must 
be supplemented to yield a determinate set, (2) this set must be solved 
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for the shock-wave boundary -value problem, thus determining n, u, and 
so forth as functions of x, and ( 3 ) the distribution function, now known 
as a function of x, must be applied to calculate the wire -temperature 
shock profile. 

The three theories used here differ, for the immediate purpose, only 
in the way they supplement the conservation equations to find a determi- 
nate set. Two basic kinetic -theory methods are involved, these being the 
Chapman-Enskog method and the Grad method. The Chapman-Enskog method 
yields first the Navier-Stokes and then the Burnett expressions for vis- 
cous stress and heat flux, as successive approximations in a series 
expansion in the mean free path. The coefficients of this series involve 
only n, u, T, and their spatial derivatives (ref. 7)* The Grad method 
relates p- i and to n, u, p, and T by taking higher order 

moments of the Boltzmann equation. This introduces further unknown vari- 
ables, and the heart of the Grad method is the choice of some physically 
sound way to truncate the procedure. The truncation which corresponds to 
the simple distribution function (b) yields the 13 -moment equations, in 
which p- • and q. appear explicitly (ref. 8). 

j-d x 

The differential equations to be solved for the shock-wave structure 
are thus the Navier-Stokes equations, Burnett equations, and 13-moment 
equations. Previous authors, notably Grad (ref. 6), Gilbarg and Paolucci 
(ref. 9)j and Zoller (ref. 10) have indicated how these equations may be 
integrated numerically, and their methods were applied quite directly. 

Once the distribution function is determined throughout the shock, 
it is used to compute the wire temperature according to the procedures of 
Bell and Schaaf (ref. 11) For the present range of shock strengths it 
is sufficient to use the approximate form of their result 

^ = Hots) + Hl (s)?S + „ 2 ( s) 53t + H 5 (s)2S ( 5 ) 

T P p J pu 


^"This method makes use of the Grad distribution function, in which 
p. i and q. enter as parameters. The Grad distribution function is 

X J X 

valid for both the 13-moment and Navier-Stokes equations. However, the 
Burnett distribution function contains higher order terms in the peculiar 
molecular velocity C, which possibly makes a significant contribution to 
the wire temperature, but to the authors f knowledge, the explicit form of 
the complete Burnett distribution function has never been worked out. 

The procedure employed in this report was to calculate Burnett wire- 
temperature profiles using the Grad distribution function, with values 
of p-i ± and q. obtained from the Burnett shock-wave solutions. No 

quantitative estimate of the error involved in this procedure is avail- 
able at the present writing. 
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for a wire lying along the z-axis with flow in the x-direction. In this 
equation, the H^(s) are functions of the molecular speed ratio, 

s = u//2RT, as given in reference 11. For a on 3 -dimensional flow of a 
monatomic gas the Stokes relation 

Pxx + p yy + p zz = ® 
is valid, as is the equation 

Pyy = p zz ~ _ 2 P xx 

regardless of the particular form employed for p when the flow is in 
the x direction. 

Finally, since radiation heat losses were obviously important in 
these experiements, the theoretical curves were corrected for these losses 
by the method of appendix A. It is interesting to observe that the radia- 
tion losses cause a decrease in maximum shock slope, as seen in the wire- 
temperature profile, but cause virtually no change in measured shock 
thickness . 

The recent kinetic theory of Ikenberry and Truesdell (refs. 12 and 13) 
has not been used in this report but is currently being applied to the 
shock-wave problem by others. Mott-Smith's theory, used in reference 1, 
was judged too unsuccessful in that work to be employed again here. 

The particular calculations which were made are listed below. 

(1) Navier-Stokes theory: T w /T q profiles were calculated for 

Mq = I. 335 , 1.454, and 1.715 by numerical integration, using the 

viscosity-temperature law given in reference 1^4. These profiles are shown 
in figures 6 (a) to 6 (d), respectively. The velocity profile at Mq = 1.576 
is shown in figure 7 (a); here w is the nondiirensional velocity defined in 
equations (Bl 8 ) . 

(2) 13-moment theory: Wire -temperature profiles were calculated for 

Maxwellian molecules at Mq = 1.335 and 1.576. These are shown in fig- 
ures 6 (a) and 6 (c), respectively. Velocity prefiles for Mq = 1.576 and 

1.713 are shown in figures 7 (a) and 7 (b), respectively. The 13-moment 
profiles become double -valued for Mq > I. 65 . 

( 3 ) Burnett theory: Wire -temperature profiles were calculated for 

Maxwellian molecules at Mq = 1.342, 1.454, 1 .^ 76 , and 1.732. These are 
shown in figures 6(a), 6(b), 6(c), and 6(d), respectively. The velocity 
profiles for Mq = 1.576 and 1.732 are shown in figures 7 (a) and 7 (b), 
respectively. 
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(4) Effect of viscosity-temperature law on theoretical profiles: 

In addition to the calculations listed above , a Navier-Stokes profile 
for Maxwellian molecules was calculated for Mp = 1.713* and a 13 -moment 
profile for the empirical viscosity law p = p(T) was calculated at 
Mp = 1,576. These profiles were compared respectively, with the Navier- 
Stokes p(T) profile at Mp = I .713 and the 13 -moment Maxwellian mole- 
cule profile at Mp = I. 576 . When these profiles are drawn with y or 
x /l* as the dimensionless distance variable, the maximum slope is 
virtually independent of the vis cosity-versus -temperature law (the point 
at which T = T* is usually very close to the point of maximum slope). 

The Maxwellian molecule profiles are slightly steeper upstream of the 
maximum-slope point, and more gradual downstream, than are the p = d(T) 
profiles. However, the differences are too slight in the present Mach 
number range to be seen in figures 6(c) and 6(d). Had x/Ap, with Ap 
the mean free path upstream of the shock, been used as abscissa, a notice- 
able difference between these profiles would be seen. Then the Maxwellian 
molecule profile would be the thickest of the family obtained from the 

viscosity laws p oc (O^ojSI), 


DISCUSSION 


Agreement Between Data and Downstream 
Rankine-Hugoniot Conditions 

It is seen from figures 6 (a) to 6 (d) that, except in the case of 
the Mp = 1.713 shock wave, the downstream end points indicated by the 
data are considerably higher than the theoretical curves. Higher values 
of T w /T 0 correspond to higher Mach numbers, and the data indicate that 
the pressure rise through the shock waves was not so great as that pre- 
dicted by theory. These discrepancies are probably due to the influence 
of the shock holder. In all the shock waves tested a region of reaccelera- 
tion followed the leveling-off region of the subsonic tail. Apparently 
the lower pressures in the accelerating subsonic flow within the shock 
holder propagate upstream into the subsonic tail of the shock wave, broad- 
ening and weakening it. The effect is less pronounced for stronger shock 
waves, as evidenced by the better agreement found here at higher Mach num- 
bers and also by comparison with Sherman's helium data. (The effect men- 
tioned above was found, to a lesser extent, also in Sherman's tests but 
it is not so apparent in the normalized profiles.) 


Comparisons With Various Theories 

Of the theoretical curves shown in figures 6(a) to 6(d), those found 
from the Navier-Stokes equations seem to fit the data best, particularly 
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in the upstream half of the shock wave. There is, however, little choice 
between this fit and that afforded by the Burnett equations; indeed the 
latter equations seem slightly more successful in the downstream half of 
the shock. This last fact is not too significant, because of the previ- 
ously mentioned possibility that the downstream tail of the shock may be 
distorted to an unknown degree by the failure to produce experimentally 
the proper downstream boundary condition. The 13-moment equations give 
a definitely inferior fit to these profiles. 

The rather close coincidence of the Navi er -Stokes and Burnett pre- 
dictions came as a surprise to the present authors, largely because a 
casual reading of Zoller's paper (ref. 10) had left the impression that 
he had already integrated the Burnett equations, finding results consid- 
erably closer to those of the 13-moment method than to those of the 
Navier-Stokes method. In fact, Zoller's basic equations do not arise 
from the Chapman-Enskog method, and while they are quite similar to the 
Burnett equations (as given in up-to-date form in ref. 15 ) they are 
evidently significantly different when applied to the shock-wave prob- 
lem. Mrs. Chang’s calculations with the correct Burnett equations 
(ref. 16 ) failed to warn of this fact because cf their poor convergence. 
On the other hand, Grad reexpanded her results and gave, in a note "added 
in proof" to reference 6, approximate formulas which predicted that the 
relative shock-wave thicknesses according to tie three theories consid- 
ered here would be just about as we have found them numerically. 

Various authors have suggested that an approximate solution by 
series-expansion techniques is most appropriate to the Navier-Stokes, 
Burnett or 13-moment equations, since these equations are themselves 
derived by series methods applied to the Boltznann equation. Good 
examples of this method are found in references l6 and 17 . 

As was done by Wang Chang (ref. l6) a series solution was con- 
structed for the shock-wave problem, for the Btrnett equations, and for 
the "third approximation" equations which contain, in addition to the 
Burnett terms, the linear third-derivative terns which arise from the 
next Chapman-Enskog approximation. The expansion technique was that 
suggested by Grad (ref. 6, p. 284) carried to crder e^, at which stage 
the first contributions of the third approximation appear. These cal- 
culations verified Grad’s formulas for the shock-wave thicknesses 
(ref. 6, p. 300) and provided the associated stock profiles. While the 
series solution of the Burnett equations, carried only to this degree 
in e, did not agree closely with the computing -machine solutions, it 
is believed that the series results predict the correct qualitative 
effect (a further thickening of the shock) of the third-derivative terms. 
The point which appears significant is that successive approximations in 
the Chapman-Enskog scheme seem to move the theoretical predictions far- 
ther from, rather than closer to, the experimental data. In this connec- 
tion it is interesting to recall that the Burnett equations can be 
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obtained, by iteration on the 13-moment equations (ref. l8, p. 258). If in 
this sense the Burnett equations are regarded as approximations to the 
13-moment equations, the fact that the 13-moment equations predict thicker 
shock waves than the Burnett equations is consistent with the aforemen- 
tioned effect of higher-order approximations In the Chapman -Enskog scheme. 

There seems to be no satisfactory explanation of the observation that 
presumably more accurate solutions of the Boltzmann equation fail to yield 
increasingly accurate predictions of the experimental data. The following 
possibilities can be suggested: (l) The Boltzmann equation is correct, 

but existing series methods of solving it are very slowly convergent or 
perhaps only asymptotic; (2) the series solutions are adequately conver- 
gent, but the Boltzmann equation itself is inapplicable in shock waves 
of moderate strength; or (3) both the Boltzmann equation and these methods 
of solution are wrong for this problem. At least, a clear indication that 
distribution functions of the form given in equation ( 4 ) are not correct 
for all shock strengths is exhibited in reference 19. There these func- 
tions are shown to become negative in certain ranges of molecular velocity 
for shock waves only slightly stronger than those investigated here. 

This is an obvious violation of the definition of a distribution function. 


Maximum-Slope Thickness 

Most investigators of shock-wave structure have used the Prandtl 
maximum— slope velocity thickness as a means for comparing the results 
of different calculations. In units of y, this velocity thickness is 
defined by 


4. u) (y) 



or, in units of x. 


4 u) w 


35 L* 

4 e 


4 u) (y) 


( 7 ) 


( 8 ) 


This is, of course, only one of several maximum- slope thicknesses that 
can be constructed. For instance, temperature or density thicknesses 
can be determined from the shock-vave solution, and from the experiments 
performed here vire-temperature thicknesses can be obtained. Most of 
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the theoretical shock-wave analyses have been presented in terms of 
velocity thickness; therefore it will be advantageous to represent the 
wire-temperature maximum-slope-thickness data in terms of equivalent 
velocity thicknesses. 

It may be remarked that a maximum-slope thickness (of any of the 
quantities, velocity, temperature, etc.) is not a particularly good meas- 
ure of the prediction of a given theory, since it depends on a purely 
local property of the shock wave. Conceivably; two different theories 
could predict similar maximum-slope thicknesses yet give profiles which 
were markedly different in shape. Both Grad (ref. 6) and Gilbarg and 
Paolucci (ref. 9) have remarked on this, and Grad introduced an "area 
thickness" which reflects better the overall shape of the profile than 
does the maximum-slope thickness. However, the area thickness is a 
rather awkward quantity to use, so the maximum- slope thickness has been 
employed here . 


The experimental wire -temperature thicknesses 


,(T W ) 


m 


exp 


obtained 


by drawing a maximum-slope line through thq dal a, were transposed to 
experimental velocity thicknesses through the correspondence 
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i 
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This ratio ran from 1 to about 1.07 as Mq increased from l.jh to 1.73; 
according to the Burnett theory, for which the most detailed profiles 
were calculated. Maximum-slope shock thicknesses transposed in this 
fashion are shown in figures 8 and 9 and tabulated in table I. In the 
figures they show ranges of ±5-percent error, which are felt to give a 
fair estimate of the accuracy with which the e>perimental maximum slope 
could be found. 


Comparison With Other Experimertal Data 

There still appears to be just one other source of experimental 
data on shock thicknesses, the work of Hornig snd his colleagues at 
Brown University (currently being continued at Princeton). This 
utilizes the optical reflectivity technique, ir which light of a wave 
length chosen nearly equal to the anticipated shock thickness is inter- 
cepted obliquely by a shock traveling down a tube. 'Given a theoretical 
density profile, and hence an index-of-refraction profile for the shock 
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wave, one can compute its reflectivity for light of a given wave length 
and angle of incidence. The theoretical result is mostly sensitive to 
the ratio of shock thickness to light wave length and comparatively 
insensitive to the detailed shape of the density profile. Hence if a 
reasonable profile shape may be assumed, the calculation may be performed 
with the shock thickness carried along as a scale factor, to be determined 
empirically by matching predicted and observed reflectivities. This work, 
including some recent and refined data on shock thicknesses in argon, has 
been summarized by Hornig in reference 20. Results shown in figures 8 
and 9 agree remarkably well with the present data, considering the com- 
pletely different experiments. Unfortunately, the data in reference 20 
attributed to Talbot and Sherman contained a small calculation error when 
supplied to Dr. Hornig. The fact that Hornig 1 s data and the present 
results agree better in figure 9 than in figure 8 is not an accident. 

The viscosity of argon is very poorly represented by p °c 
(a) = Constant) over any considerable temperature range. In the present 
investigation wind-tunnel data T* ~ Constant = 225° K, and 53 0.9, 
while in Hornig' s data, T]_ » Constant = 300° K, and o> » 0.8. Fig- 
ure 8 gives the truer comparison of these two sets of data since &m/l>* 
is virtually independent of moderate change's in the viscosity- temperature 
law. 


CONCLUDING REMARKS 


The profiles and thicknesses of normal shock waves in argon at Mach 
numbers of 1 . 335 , 1 - 576 , and I .713 were determined experimentally 

by means of a free-molecule probe whose equilibrium temperature is 
related by kinetic theory to the local flow properties and their 
gradients . 

Comparisons were made between the experimental shock profiles and 
the theoretical profiles calculated from the Navier-Stokes equations, 
the Grad 13-moment equations, and the Burnett equations. New, very 
accurate, numerical integrations of the Burnett equations were obtained 
for this purpose, with results quite different from those found by 
Zoller, to whom the solution of this problem is frequently attributed. 

The experimental shock profiles were predicted with approximately 
equal success by the Navier-Stokes and Burnett theories, while the 
13-moment method was definitely less satisfactory. A surprising feature 
of the theoretical results is the relatively small difference in pre- 
dictions from Navier-Stokes and Burnett theories in the present range 
of shock strengths and the contrastingly large difference between pre- 
dictions of the Burnett and Grad theories. It is concluded that the 
Navier-Stokes equations are correct for weak shocks and that within the 



20 


present shock-strength range the Burnett equations make no improvement 
which justifies the trouble of solving them, lor shocks of substantially 
greater strength, it remains fundamentally doubtful that any of these 
theories can be correct. 


University of California, 

Berkeley, Calif., January 1, 1958* 
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APPENDIX A 

FREE -MOLECULE -PROBE HEAT LOSSES 


Radiation-Convection Heat Balance 

Consider first a probe in a uniform stream. The probe has an equi- 
librium temperature T w which is less than the adiabatic equilibrium 
temperature T aw . Heat is transferred to the probe by convection and is 
lost from the probe by radiation. (For the time being, conduction losses 
and i^R heating will be neglected.) 

A heat balance for the probe is, per unit length, 

2h c itra(T av - T w ) = 2 itra^e w T w ^ - CqTo^^ (Al) 

in which e w is the emissivity of the wire at temperature T w and e 0 
is its emissivity at the surrounding temperature T 0 . For a pure metal, 
e w/ e o = VT 0 (ref. 21 ), so 


h c^(T aw 


- T w ) = 


<^oV 



(A2) 


If this equation is used to correlate the free-stream equilibrium- 
temperature data, with 


lie 



(1 + s2)l 0 



+ 



(A3) 


then it is found that the single value e Q /a = 0.13 gives a curve for 
T v /T 0 which passes through all of the data points, as shown in figure 5 . 
In this figure, the Taw/^o curve was calculated from the functions 
given in reference 5 * 

The value e Q / a =0.13 seems to be rather large. A typical value 
for the emissive power of tungsten at room temperature is e G = 0 . 03 , 
and values cited in reference 22 for the accommodation coefficient of 
argon on tungsten range from a = 0.55 (clean tungsten) to a = 1.0 
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(surface saturated with argon). However, it is not clear whether the 
above values apply for the situation at hand. For example, there was 
undoubtedly a layer of adsorbed air on the wire, and also the surface 
condition of the wire, which affects both e Q and a, was unknown. In 
addition, the wire geometry itself was rather complicated, since it was 
constructed with two copper-plated lap joints at the potential lead junc- 
tions, and these double-strand joints make up almost 20 percent of the 
measuring span. Because of these aforementioned uncertainties, it seems 
most reasonable to accept the value e 0 /a = 0.1;> simply as a calibration 
constant and to refrain from more specific interpretation. 


Estimate of Conduction End Losses 

The full heat balance for the wire, including conduction and current 
heating but neglecting changes in resistance anc. conductivity with tem- 
perature, is 


d%\ 


w 


dx 


2ae 0 1'o 

Kr 


SeT - 1 

.tJ 


2h c a / . i_2p> T 

+ (^aw " ^w) + ~ 0 


Kr 


(A4) 


jtKr 


The boundary conditions at x = 0 are, wiih x = 0 taken at the 
center of the wire span, 


dL 

dx 


w 


0 


and, at x = ±Z, 


T 


= T 


o 


Putting 


T w “ T 0 (l + t) 


^aw “ + *^aw) 
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equation (A4) may be linearized and solved, with the result 


t 



cosh cx 
cosh cZ 


(A5) 


where 


b 


2h c at av 

Kr 


i 2 R 

*Kr 2 T 


I 


o 


5 

10ae o To 2h c a 

c = + 

Kr Kr 

For a numerical example , consider "the conditions which were obtained 
at M]_ = 1.713 in the free stream. Take 

a = 5.6733 X 10 "5 ergs/cm^ °K sec 
K = 2 X 10T ergs/cm °K sec 

r = 3.17 X lO - ^ cm 
T 0 = 296.5° K 

i = 1 X 10”^ amp 
R 1 « 15 ohms / cm 
OUp ~ 0.00332/°K 

h c = 1.51 X 10^ ergs/cm 2 °K sec 
e Q /a = 0.13 
a = 1.0 

^aw/^o = 1-221 
*^aw “ 0.221 
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Then, c =2.25 cm"! and — = 0.2075* It is found also that the effects 

of i^R 1 heating are negligible. The most difficult thing to estimate 
is the effective half length Z between the wire center and the assumed 
end-point heat sink at temperature T 0 . In the actual tests the current 
and potential leads spanned about 19 centimeters . However, only about 
10 centimeters of this span was immersed in the uniform stream, the 


remaining 4^ centimeters on each end being in ;he nozzle boundary layer 


and outside the jet. It is known also that th 2 equilibrium temperature 
of the wire within the boundary layer is actually higher than it is in 
the uniform flow, so there is in fact some hea: input to the wire in 
these regions. Since only an estimate of the conduction effect is 
needed, the rather pessimistic choice cZ = 10 is made. Then, 


t 


(x 


0.2075(1 


cosh 2.25-X 

11,000 " 


(A6) 


The measuring span of the wire was about l centimeters, so with 
x = 1, 


(t) x=1 » 0.2070 (A7) 

Even if the conduction losses are doubled to account for heat flow out 
the potential leads, it is evident that these Losses are still negligible, 
since the difference between the center-point temperature and the tempera- 
ture at x = 1 is, from equation (A7), 

(Tj x=0 - (T w ) x=1 - 0.00l(T aw - T 0 ) - 0.065° C 

Note that the center point temperature ratio with correction for 
radiation obtained from this linearized calculation is 


= 1.2075 

x=0 

which is quite close to the exact value obtained from equation (A2) 
for the flow conditions chosen. 
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Radiation Corrections to Theoretical Shock Profile 

For comparison with the experimental results, the calculated ^av/^o 

profiles were corrected, point by point, for radiative heat losses. The 
heat-transfer coefficient which is appropriate to a nonuniform flow was 
taken from reference 11 and is 



This heat-transfer coefficient and the value e Q /a = 0,13 
with equation (A2) to obtain the corrected 'profiles . 


were employed 
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APPENDIX B 

INTEGRATION OF BURNETT EQUATIONS 


The Burnett equations for one -dimensional flow are nearly identical 
to those integrated for a shock wave by Zoller in 1951* Zoller's method 
of integration is so thoroughly presented in reference 10 that it was 
quite simple to set up the Burnett equations fcr the solutions described 
here by merely copying his procedures. This wes done, as follows. 

The three conservation equations for one-dimensional steady flow of 
a fluid are 


j-(pu) = 0 
dx 


\ 



+ P + Pxx) 


= 0 


y (bi) 


a_ 

dx 


pu 



+ u(p + Pxx) + c x 


= 0 


/ 


in which p is the mass density, u, the flow velocity, p, the pres- 
sure, E, the internal energy, q x , the x component of heat flux, and 
Pxx> the x component of excess viscous pressure. For a perfect mon- 
atomic gas, 


E = | RT 
p = pRT [ 


(B2) 


The differential system of equations (Bl) may te integrated once to yield 


\ 


pu = m 

pu 2 + p + Pxx = p > 

pu3 + 5pu + 2p xx u + 2q x = (1 


(B3) 


where P, Q, and m are constants of the flow 
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The viscous stress may be expressed as 

— 1 , tt 

PxX - 1*XX + Pxx 


(B4) 


with p ' = - — (j, — , the familiar Navier-Stokes expression. 


4 

- d 

3 dx 

the heat flux q x may be written 


q x = q x' + q x" 


Similarly, 


(B5) 


dT 

with q ' = -k — . the Fourier heat- conduction law. 
x dx 


Wang Chang and Uhlenbeck (cited in ref. 15) give the currently 
accepted expressions for the Burnett stresses and heat flux, which reduce 
for one -dimensional flow to the following equations: 


■fa 


Ik 
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+ OV 

9 d 27 ^ 


+ 2 E 2 ^ + S 
3P T ^dx 2 3 




H_ <*E dT + 
ppT dx dx 


2 ^ )g_(*\ _ 2 £ ^ dfp + 2 

3 


dp dp 

Ok 1 — \ - — ’■ OUo — + — 6L>) — 

3 pp ^ ^.,2 3 dx dx 


pT 2\dx/ 3 PP dx 2 5 p 2p 

q/ = idi + 8 e 2 + 2 9 2 e W ^ + 2 p£/ ^ _ 

x ptI 1 3 3 5 3 7dxdx 3 p\ ^dx 2 


(B6) 


2 0 , 

3 ? pu\dx y 


The values of the constants 9^ and 

= 10/3 
0^2 = 2 
0*5 = 3 
= 0 
= 3 

“6 = 8 


o>i are, for Maxwellian molecules 
01 = 75/8 
e 2 = -45/8 
e 3 = -5 

04 = 5 

0 5 = 117 A 
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These numbers are now introduced, and the continuity equation and 
equation of state are used to eliminate the pressure and density vari- 
ables. This gives the following set of six equations for six unknowns: 


. mRT , , tt „ 

mu + + p + p vv = P 

u ^xx ^xx 


mu 


2 + 5mRT + + p xx") + 2(q x ' + q x " ) = Q 


Pxx 


it n 

3 dx 


q 1 = -k AS 

dx 


> (B7) 


■ XX 


v 2 2 

40 P^ / du \ _ 4 p£/du\ + 4 p£ dfu { 4 p 2 du dT + 

27 mRT \dx / 3 mu\dxy 3 m ^2 3 mT dx dx 


2 p£ u dfT 

3 m T ^,,2 
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dx mi 


nm 

ip2 ^dx 




Ir 95 P^ u du dT 7 P^ d^u 0 p2/du\ 

q = — - - — u + 2 — — 

x 8 m T dx dx 4m ^2 m \dx f 
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The dimensionless variables suggested by Zoller are 


_ mj/2RT 


Pxx' = u Pxx' 
p mRT 


Pxx = u Pxx 
p mRT 


II U 
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4 lx u_ 

5 mRT 


3 p m 

z — — t^- — x 
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In these variables, the system of equatior s (B7) becomes 
2v^ - 2v + T^(l + T) t + t] m ) =0 


v 2 + t 2 


| + tT + V' - f -(!' + I") 

2 4 v 


= a 


V = 




dv 

dz 


> (B9) 


I' =2 v St 

4 dz 


v . - f| .2i£\ (v) 2 t iit( S ,)2 ix,. { . 

8 dz 6 dz y6 4 V 2J 27 6 v 


t- = - ®-Tv2H-SIl( n .)2 + w,. { . 

l 60 dz 32 v 8 


The fourth of these equations involves the assumption that the Prandtl 
number is 2/3, which is satisfied exactly by Mftxwell molecules and is 
nearly true for all pure monatomic gases. 


If in equations (B9) the first two members', are used to eliminate 
V and from the last two members, a little rearranging of terms 

then leads to the following set of four first-order, nonlinear, ordinary 
differential equations 
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dv 

_ nl 


dz 
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dT 

it 11 


dz 

9 V 


dV 
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dz 
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R x ■ 

TV 

40 

dz 

7V 2 


(BIO) 
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in which 


Hn = 


2v - 2v £ 


1 v2 6 


- 1 - v + (r — - ^l(v) 2 + 7 - t i' r -^d') 2 


27 


> (Bll) 


„ 2 - i^[3( 

5t^ 


t 2 - t v 2 ) + 4(v - v v ) - 2(v 2 - v v 2 ) - 6' - ^ i'V + ^ IC 1 !') 2 


The quantities r v and v v , which appear in Up, axe the boundary- 

values of t and v far upstream or far downstream of the shock, that 
is, at z = +bo. At these locations the stress and heat flux are pre- 
sumed to vanish, so that the first two members of equation (B9) can be 
quickly solved for t v and v v in terms of a, the shock strength 
parameter. This gives first 


2v v 2 - 2v v + t v 2 




(B12) 


The index v will be given the value 0 to indicate conditions far 
upstream of the shock (z = -») , and v = 1 to indicate downstream 
conditions . 


Since the boundary conditions are specified at the singular points 
of equations (BIO), a special investigation of the solutions in the 
vicinity of these points is a prerequisite for numerical integration. 
For this purpose one sets 


v = v v + ft 

T = T v - d 

v = v 

i' = s’ 


(B13) 


where ft, i3, V , and- are quantities negligible to any order higher 

than the first in the vicinity of the singular point. Next one assumes 
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2 = Qe 




,Az 




T| ' = T]e^ z 


6' = Se' 


(Bl4) 


Then there is obtained, from equations (BIO) the set of linear algebraic 
equations 


AO = - — 


A$ 
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9 v v 

160 ( k n_ 
63t v v v U t v 


12 v v - 

— — ■d - i 

5 T 2 
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M - > - -h&v - Do - nl - -=^2 - - « - s 

V V v\ n V Ty 


40 Ik SL 12 v v 1 7 

7v v 2 l5 x v 5 T 2 


which has nonzero solutions only if A satisfies the characteristic 
equation 

l89v v ^(l - v v )^A - 1080v v ^(l - v v )5a 5 + 2328v v ^(l - V V ) 2 A 2 - 


96 v v (l - v v )( 38 v v - 15 )A + 256(8v v - 5) = 0 


(Bl6) 


From the four roots of this quartic, one which neets the requirements of 
physical reality can be selected for each value of v. The requirements 
are that the real part of A be positive for v = 0(z = -oo) and nega- 
tive for v = l(z = oo) and that A -> 0 as v o “* v q ->5/8 (that is, as 

the shock strength vanishes). 
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When A is a root of equation (Bl6) the first three members of 
equations (B 15 ) can be solved for the ratios $/R, T)/R, and |/R* 

yielding 



The initial values of v, t, tj, and | are then established for the 
purposes of computation by choosing the point z = 0 to correspond to 
some small but finite value of R. The initial derivatives at that point 
can then be found from equations (BIO) or by differentiating equa- 
tions (Bl 4 ) . The two methods will give nearly the same numbers if the 
starting value of R is sufficiently small. 

To insure the stability of machine integrations, one must determine 
by examination of tbje roots of equation (Bl6) whether to start near the 
upstream limit of the shock and integrate downstream or vice versa. The 
real parts of these roots are plotted versus v v in figure 10. 

At v = 1 (0.25 = v y ^ O.625) there are either four real roots, 

three positive and one negative, or two real roots of opposite sign plus 
a pair of complex conjugates with positive real part. The physically 
acceptable root is the negative real one. This singular point is of a 
saddle-point type and is ideal as a starting point for numerical 
integration. 

At v = 0 (0.625 v v ^ l) the singular point has a variety of 

types, none of which are suitable for the initiation of numerical integra- 
tion. These types dictate qualitatively the behavior of the numerical 
solution as it approaches its upstream limit. Thus, for 0.625 % v v < 0.80 

the solution will come in monotonically, for 0.80 < v v < 0.9^ a damped 
oscillatory approach can be expected, and for 0.9^- « v v ^ 1 this oscil- 
lation may amplify rather than damp out. This behavior has been qualita- 
tively confirmed by the computing machine, which also indicated that up 
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to Mach 2 the oscillations encountered are so small as to be barely 
noticeable. 

Integrations were consequently started from downstream, with the 
arbitrary choice 0 , = 0.001, and proceeded toward increasingly negative 
values of z. They were performed on the NAREC digital computer of the 
Naval Research Laboratory, Washington, D. C. Numerical experiments 
confirmed that, with a Runge-Kutta third-order integration method 
(ref. 23 ) and integration steps equal to about 0.025 times the maximum- 
slope shock thickness, the results were independent of further reduction 
of step size to an accuracy considerably greater than that required here. 
The results also showed a high degree of independence of small variations 
in the initial derivative data, as would be expected of an integration 
starting out of a saddle point. These computations were programed and 
supervised by Mr. John L. Hammersmith of NRL, to whom the authors are 
greatly indebted. 

To express these results in the nomenclature used elsewhere in this 
report, the molecular speed ratio s = v/r and the following relations 
were necessary: 


Pxx = 2v - 2v 2 
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^x _ a - 2v + v^ 3 
pu r 2 2 
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^25 - 16a 
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TABLE I.- SUMMARY OF FLOW PROPERTIES 
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Equivalent velocity thicknesses, obtained from measured wire -temperature thick- 
nesses by use of equation (9)> 
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Figure 2.- Shock holder and wire-traverse mechanism. 
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Figure 3.- Glow photographs of shock-holder operation. 
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Figure 4 .- Schematic diagram of resistance-wire probe. 
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Figure 6.- Temperature profiles. 
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Figure 6.- Continued. 
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Figure 6.- Continued. 
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Figure 6.- 
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Figure 7* “ Velocity profiles. 




Figure 7-~ Concluded. 
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Figure 8.- Shock-wave thickness 
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Figure 9*~ Shock-wave thickness in terms of apstream mean free path. 
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Figure 10.- Roots of equation (Bl6) . 
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